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ABSTRACT 

We present a near-infrared study of the candidate star cluster Mercer 81, located at the centre 
of the G338.4+0.1 Hll region, and close to the TeV gamma-ray source HESS 1640-465. Us- 
ing HST/NICMOS imaging and VLT/ISAAC spectroscopy we have detected a compact and 
highly extincted cluster of stars, though the bright stars in the centre of the field are in fact 
foreground objects. The cluster contains nine stars with strong Pa emission, one of which 
we identify as a Wolf-Rayet (WR) star, as well as an A-type supergiant. The line-of-sight ex- 
tinction is very large, Ay ~ 45, illustrating the challenges of locating young star clusters in 
the Galactic Plane. From a quantitative analysis of the WR star we argue for a cluster age of 
3.7+Q g Myr, and, assuming that all emission-line stars are WRs, a cluster mass of p> 10 4 M©. 
A kinematic analysis of the cluster's surrounding Hll-region shows that the cluster is located 
in the Galactic disk at a distance of 1 1 ±2 kpc. This places the cluster close to where the far 
end of the Bar intersects the Norma spiral arm. This cluster, as well as the nearby cluster 
[DBS2003] 179, represent the first detections of active star cluster formation at this side of the 
Bar, in contrast to the near side which is well known to have recently undergone a ~ 10 6 M© 
starburst episode. 

Key words: (Galaxy:) open clusters and associations: general (Galaxy:) open clusters and 
associations: individual: Mercer 81 stars: Wolf-Rayet (ISM:) H ii regions ISM: clouds 



1 INTRODUCTION 

Young massive star clusters (YMCs - ages < 50Myr, masses > 
10 4 Mq) have relevance to many areas of astrophysics. They con- 
tain large numbers of massive stars, whose high temperatures, high 
luminosities, dense winds, and supernova explosions make YMCs 
a considerable source of mechanical energy, ionizing radiation and 
chemically processed ejecta. The effect that they have on their sur- 
roundings is profound, clearing away the remains of their natal 
molecular cloud, whilst reveali ng and triggering subsequent gen- 
erations of s t ar formation (e.g. iGonzalez Delgado & Perezll2000l : 
ISmithl 120061 : iDavies et all 12011 at iDe Marchi et al.1 1201 lh . Their 
large populations of massive stars make them ideal natural labo- 
ratories with which to stud y the evolution of massive stars up to 
supernova and beyond (e .g. Martins et al. 2007, 2008; Bibb v et all 
l2008l:lDavies et al.ll2008ll2009h . FinalhTthev can dominate the ra- 
diative output of their host galaxies, either through direct ultraviolet 
and optical emission, or thro ugh reprocessed emission in the form 
of ionised gas or heated dust (Al onso-Herrero et al . 2002). 



Our knowledge of our own Galaxy's population of YMCs is 
extremely incomplete, in stark contrast to exter nal galaxies (e.g. 
iBastian et aT]|2005l : iKonstantopoulos et al.ll2009l) . The high levels 
of interstellar extinction in the plane of the Galaxy have meant that 
until recently very few clusters were known beyond a distance of 
~lkpc. Most known clusters beyond this distance w ere found by 
targeted searches of, fo r example, the Gala ctic Centre dCotera et all 
1996l : lFigeretal.l 19991) . giant Hn-regions dBlum et al.1l 1999L|2000L 



2001), fields aro und newly-born neutron stars dFuchs et alJI 199S 



Vrba et alE oOO). or simply because the foreground ex tinction was 



low e nough to detect the cluster at optical wavelengths iWe sterhindl 
Il987h . However, recent i nfrared (IR) surveys of the Galactic plane, 



beginning with 2MASS (Skrutskie et al 



Spitz er/GLIMPSE (Benja min et alj|2003l) and VVV jMinniti et all 



2006), and more recently 



2010), are at last helping us to uncover the cluster population of the 
Galactic disk, and affording the opportunity to search the Galactic 
Plane for clusters in a systematic way. 

By-eye and a lg orithmic se a rches o f survey da ta (e.g . 
Ilvanov et al.l 120021 : iDutra et al.l 120031 : iMercer et all 120051; 
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iFroebrich et al.l 120071 : iBorissova et alj 1201 lh have yielded over 
1,000 candidates for newly discovered star clusters. However, such 
catalogues inevitably contain large numbers of false positives, due 
to chance alignments of stars, patchy foreground extinction, and 
spatially extended emission incorrectly classified as unresolved 
star clusters. Therefore, these catalogues of candidates must be 
analysed carefully using multiwavelength survey data and, ulti- 
mately, follow-up spectroscopic observations b efore their distances 
and physical properti e s may be derived ( e.g. iKurtev et alj|2007l ; 
iMessineo et al.l 120091 : lHanson et ail l20ld) . Only then can they 
be placed in the framework of the Galaxy's recent star-forming 
history. 

One such c andidate is object #81 from the catalogue of 
iMercer et alj d2005t) . known hereafter as Mc8lQ. This object was 
found in an algorithmic search of the GLIMPSE survey, by look- 
ing for spatial groups of stars with similar photometric properties. 
The clues to its nature, however, come from cross-correlation with 
other data in the literature. The object appears to be at the centre 
of the Hll-region G338. 4+0.1, a bubb le of warm dust and ionized 
gas, visible in the SU MSS 843MHz dBock et al.|[l999h and MSX 
8/xm jPrice et al.l2 001) images, and in more detail in the GLIMPSE 
8/im image (Fig. [T V The o bject' s location is close to a supernova 
remnant (SNR) by iGreer] d2004l) from the 'shell' morphology of 
continuum radio emission, but with no spectral index measurement 
the object could also be a wind-blown bubble. 

Close to the centre of the S NR is a high-energy Te V gamma- 
ray source, HESS J 1640-465 dAharonian et alj l2005h . at a dis- 
tance of 8' (a linear distance of 22pc, if the complex is at a dis- 
tance of llkpc - see Sect. |3.3t . This source, as with many other 
TeV sources, is thoug ht to be associated with a pulsar wind neb- 
ula dFunk et aT 2007), indicative of recent SN activity. Indeed, 
TeV emission can be a useful tracer of massive star formation 
in the Galactic Plane as it is unaffected by absorption, and there 
are other young massive star clusters in the literatu re that are 
known to be associated with such sou rces - RSGC1 dFiger et al.1 
l2006l:lDavies et alj2008h, CI 1813-17 8 dMessineo et alj2008h . and 
Westerlund 2 dAharonian et alj|2007l) . Finally, in a follow-up of 
HESS J1640-465, a number of hard X-ray sources were detected 
in the field of Mc81, with on e source being perfectly aligned with 
the cluster dLandi et al . 2006, see right panel of Fig.[JJl. Such emis- 
sion could be explained by either a recently-formed neutron star or 
a colliding wind binary, with both explanations being indicators of 
youth (for an analysis of the X-ray emission f rom another young 
star cluster, Westerlund 1, see lClark et~ai]|2008h . 

Based on this evidence we have extensively followed-up Mc8 1 
with near-infrared (NIR) photometry and spectroscopy, with a view 
to confirming that the object is indeed a young star cluster, and 
ultimately to determine the cluster's physical properties. 

We begin in Sect.|2]with a description of the observations, data 
reduction and analysis steps. In Sect.[3]we describe our results, and 
show that the object is indeed a highly extincted star cluster, and 
estimate its age and mass. We summarize our results in Sect. [5] 



1 Objects in this catalogue are referred to by some authors with the prefix 
GLIMPSE, e.g. GLIMPSE81 



Table 1. Read-sequences and total integration times employed for each filter 
during the NICMOS observations. 



Filter 


SAMP-SEQ 


NSAMP 


T int (S) 


F160W 


STEP2 


15 


144 


F222M 


STEP8 


12 


336 


F187N 


STEP8 


10 


240 


F190N 


STEP8 


10 


240 



2 OBSERVATIONS & DATA REDUCTION 
2.1 Imaging 

Images were obtained with HST/NICMOS on 22 October 2008, 
as part of observing programme #11545 (PI: B. Davies). We used 
the NIC3 camera which has a field-of-view of 5 1. 2" x 5 1.2" and a 
pixel scale of 0.2". We observed the cluster through filters F160W 
and F222M, as well as the narrow-band filters F187N and F190N 
which are centred on Pa and the neighbouring continuum respec- 
tively. In addition to the cluster we observed a nearby control field 
through the F160W and F222M filters in order to characterize the 
foreground population. The fields of observation are indicated in 
Fig.E 

Our observations used a spiral dither pattern with six pointings 
with offset distance was set to 5.07". This sub-pixel dithering tech- 
nique was designed to minimise the impact of non-uniform intra- 
pixel sensitivity on our photometry. The MULTIACCUM read 
modes were used, with the sampling sequences and total integra- 
tion times that are listed in Table[JJ 

Our reduction procedure followed the guidelines of the NIC- 
MOS Data Handbook v7.0. The standard reduction steps of bias 
subtraction, dark-current correction and flat-fielding were per- 
formed using CALNICA. Before mosaicing, each dithered observa- 
tion was subsampled onto a 3 x finer grid using bi-linear intepola- 
tion to account for the sub-pixel dithering. 

Photometry was extracted from the images using th e 
STARFINDER package which run within IDL dDiolaiti et"all 2000). 
in conjunction with point-spread functions (PSFs) which were 
computed for each filter using TINYTIM. STARFINDER uses these 
PSFs to locate stars within each image, and we employed two iter- 
ation cycles to fine-tune the astrometry and photometry. Since our 
fields of observation are not very crowded, we did not use the de- 
blending algorithm, as this was found to produce many false detec- 
tions. Uncertainties and completeness limits were determined by 
inserting fake stars into each image and measuring the recovery 
rate. 



2.2 Spectroscopy 

Spectroscopic data were taken on the nights of 2009-4-11 and 
2009-5-4 as part of ESO observing programme 083.D-0765(A) 
(PI: E. Puga), using the ISAAC spectrograph on the VLT 
dMoorwood etaill 19981) in 'medium' resolution mode. Our targets 
were the stars labelled '1' and '3' in Fig. [2] which were deter- 
mined to be likely cluster members based on their photometric 
properties (see later). We used the 0.8" slit at three different cen- 
tral wavlengths: 1.71/im, 2.09/im, and 2.21/im, providing a spec- 
tral resolution of AA/A ~4,000. The DITxNDITxNINT combi- 
nation for each wavelength setting was (8x8x30s), (8x8x32s), 
(16x8x24s) respectively. The observations were taken in a ABBA 
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MSX8[im IRACS^m 




Gal I Offset (arcmin) Gal I Offset (arcmin) 

Figure 1. Wide-field images of Mc81 and its surround ings. Left: the MSX 8/tm image, which shows the diffuse nebula surrounding the cluster. The locations 
of Swift X-ray sources are indicated (Landi et al. 2006), as well as the location of the TeV emission source HESS J1640-465. The contours indicate the 
morphology of the 843MHz emission IBock et al.ll 199911 . and the dashed green box shows the field-of-view of the right-hand panel. Right: a higher resolution 
8fj,m image from the GLIMPSE survey. The brightest stars of the cluster are coincident with the X-ray source Sw 2. The NICMOS cluster and control fields 
of view are illustrated by the dashed green boxes. 
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Figure 2. NICMOS images of the cluster. Left: image of the cluster taken through the F222M filter. The two stars for which we have spectra, as well as the 
other emission-line stars, are indicated by the blue triangles. Right: the difference image (F187N-F190N), which highlights the emission-line stars. The arrows 
in the top-right of each image indicate the orientation. 



pattern to isolate and subtract sky emission features. In addition 
to the target stars we observed the B9v stars Hip090248 and 
Hip091286 as measures of the telluric absorption, as well as the 
usual observations of flat-fields, dark frames and arcs for wave- 
length calibration. To characterize the spatial distortion, a bright 
field star was stepped along the slit and re-observed multiple times. 

The data reduction procedure began with subtraction of nod 
pairs to remove sky emission, and division by a normalized flat- 



field. The 2-D frames were then rectified onto an orthogonal grid, 
using the stepped-star and arc frames to characterize the dis- 
tortion in the spatial and dispersion directions respectively. This 
process also wavelength calibrates the data, with r.m.s. residuals 
which were found to be less than a tenth of a resolution element 
(~10km s _1 ). After rectification, the spectra were extracted and 
combined. 

The telluric standard spectra had their intrinsic H I absoption 
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Figure 3. Three-colour image of the cluster from the NICMOS data, with colours as follows - R=F222M, G=F160W, B=(F187N-F190N). Therefore, the 
highly reddened cluster stars appear red/orange, and the cluster's emission-line stars appear pink/megenta. 



removed by fitting the lines with Voigt profiles. The telluric spectra 
were then cross-correlated with the target spectra in the region of 
isolated telluric features to correct for any residual sub-pixel shifts, 
before the target spectra were divided through by the telluric spec- 
tra. 



3 RESULTS & ANALYSIS 
3.1 Photometry 

The NICMOS images of Mc8 1 can be seen in Fig.[2] The left panel 
shows the F222M image of the cluster. In the right panel, we show 
the difference image [F187N-F190N], which clearly highlights 9 
stars with significant Pa emissioifl This strongly suggests that 
these are hot stars with strong winds, and are therefore likely to 
have ages <Cl0Myr. 

In Fig. [3] we show a 3-colour RGB image of F222M (red), 

2 It is possible that He II emission contributes to the flux in this band, es- 
pecially if the stars are WRs. 



F160W (green), and [F187N-F190N] (blue). Around the coordi- 
nate centre, there is a clear group of highly reddened stars which 
form the putative cluster. In this colour scheme, the emission line 
stars, which are also heavily reddened, show up as pink/magenta. 
From these data there is already strong evidence for a young, highly 
reddened cluster of stars in the field of Mc81. The four bright stars 
to the south of the cluster have a green/yellow colour, indicating 
that they are in the foreground. Ironically, it is likely that these 
four stars in part tri ggered the cluster detection algorithm used by 
iMercer et alj d2005l) . These authors claim that they detected an as- 
sociation of 65 stars, though the positions of these stars are not 
listed, so it is not possible for us to check whether the cluster de- 
tected by Mercer et al.'s algorithm bears any relation to the cluster 
we describe here. For the sake of clarity and consistency, we con- 
tinue to refer to the cluster studied in this work as Mc81. 

The results of the NICMOS photometry are shown in Fig. 
[4] The left panel shows a colour-magnitude diagram of the stars 
within 18" of the cluster centre, which we define as the location 
of the brightest emission-line star. Also shown on the plot are data 
from an area in the control field of identical size. There is an ob- 
vious difference between the two regions, with the 'cluster' field 
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Figure 4. Colour-magnitude diagrams for Mc81 from the NICMOS photometry.Le/f: The stars within 15" of the cluster centre, which we define as the position 
of Mc81-2, compared to stars in a nearby control field of the same angular size. Centre: the Pa emission of stars within 15" of the cluster centre, as determined 
from the colour (mig7 — migo). Stars with significant emission are marked with green circles. Right: the same as the left panel, after the cluster field has 
been decontaminated of foreground stars using the control field observations. The long-dashed lines in each figure show the 50% completeness levels. 



showing an excess of stars at (mi6o — 771222)— 2.3. The centre panel 
shows the Pa excess stars. Those stars with (mi87 ™ w.190) colours 
<0.3 and 771222 >13 are defined as 'Pa emitters', and are indicated 
on the plot. The locations and photometry of these stars are listed 
in Table|2] 

Finally, the right-hand panel of Fig. [4] shows the same as the 
left, after the cluster field has been decontaminated of field stars 
using the data from the control-field. We do this by eliminating 
stars from the cluster field which have a control-field star close 
by in colour-magnitude space. We set this limiting distance to be 
the larger of either the cluster field star's 1-a photometric errors, 
or 0.15mags in colour and O.lmags in magnitude. We also set the 
restriction that no control-field star can eliminate more than one 
cluster field star. After decontamination, the cluster sequence at 
(mieo — 777,222)— 2.3 becomes clearer, though there is still some 
scatter. 



Table 2. Astrometry and photometry of Mc81-1, plus the emission-line 
stars. Astrometry is taken from the HST observations, and comparisons to 
2MASS indicate that it is accurate to 



ID 


RA-DEC (J2000) 


"1160 


777222 


mis7 


"7190 


1 


16 40 29.83 


-46 23 33.9 


11.14 


8.90 


9.67 


9.75 


2 


16 40 29.65 


-46 23 29.1 


13.25 


10.59 


10.69 


11.63 


3 


16 40 30.08 


-46 23 11.4 


12.97 


10.82 


10.75 


11.60 


4 


16 40 29.65 


-46 23 28.7 


13.44 


10.95 


11.15 


11.91 


5 


16 40 28.35 


-46 23 25.6 


14.13 


11.42 


11.67 


12.45 


6 


16 40 29.32 


-46 23 11.6 


13.46 


11.42 


11.56 


12.14 


7 


16 40 29.60 


-46 23 25.6 


13.55 


11.45 


11.56 


12.31 


8 


16 40 28.94 


-46 23 27.1 


14.27 


11.46 


11.95 


12.65 


9 


16 40 29.32 


-46 23 38.2 


13.76 


11.82 


12.16 


12.54 


10 


16 40 30.05 


-46 23 24.3 


15.08 


12.78 


13.41 


13.70 



3.2 Spectroscopy 

We now present the results of the spectroscopic obse rvations. To 
classify the spectra , we re fer to the s pectral atlases of j Hanson et al.l 
d2005l) . lFiger et all d 19971) . as well as lCrowther et all d2006h . 

The spectra of the two brightest stars in Mc81 are shown in 
Fig. [5] The brighter star, Mc81-1, has relatively weak spectral fea- 
tures. The absorption lines of the Hydrogen Brackett series are 
seen, as well as some faint Mg II emission. We see no evidence of 
He I in either absorption or emission. This absence of He I suggests 
a spectral type later than ~B5, while the emission lines of Mg II are 
indicative of a substantial stellar wind, more typical of supergiants. 
The 'ripples' seen between 2.05-2.08/iim are due to poor cancella- 
tion of the CO2 absorption. For now, we assign a loose spectral type 
to this star of late-B/early-A supergiant. Taking the distance and ex- 



tinction deri ved in Sectionl3.3l as well as the bolometric corrections 
tabulated bv lBlum et alj d2000h for spectral types B7-A2, we esti- 
mate the luminosity of the star to be in the range log(L / Lq) =5.4- 
5.8, placing the star close to the empirical stellar lu minosity limit 
atlog(L/L ) ~5.9 (Hump hreys & Davidson! 19791) . Such stars are 
also seen in Westerlund 1 dClark et al.ll2005h . The star's proximity 
to the Humphreys-Davidson limit suggests that the star may be in 
an unstable phase of evolution, such as a Luminous Blue Variable 
or Yellow Hypergiant phase, though further spectroscopic and pho- 
tometric monitoring would be required to verify this. 

In contrast to Mc81-1, Mc81-3 has a spectrum rich in strong, 
broad emission lines. These lines can be attributed to Hi, He I, 
Hell and Nm. The ratio of the Hell 2.189/j.m to the complex at 
2.115/^m, as well as the absorption of Hel 2.189/j.m, allow us to 
tightly constrain the spectral type of this star to be WN7-8. 
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Figure 5. Spectra of the two stars observed. The wiggles seen at ~2.05/-tm are due to fringing on the detector, which we were unable to remove. 



3.3 Extinction and distance 

To calculate the extinction, we first determine the reddening of 
the cluster sequence from the right-hand panel of Fig. [4] The av- 
erage colour of the stars in the decontaminated cluster field is 
(mi6o — rri222)=2.3±0.3. If we make the approximation that all 
main-sequence stars that we detect should have colours of approx- 
imately zero, the observed average colour is due to extinction. We 
can then determine the extinction from the following relation, 



Ay 



(Ai/A 2 ) c 



(1) 



with Ai —1.60fim and A2 =2.22/im, i.e. the wavelengths of 
the NICMOS F160W and F222M filters. The parameter a has 
been studied by nu merous authors in recent years (see e.g. 
IStead & Hoarefe o09. and references therein), with the most con- 
temporary measurements converging on a = —2.0 ± 0.1. This 
therefore implies that the extinction towards Mc81 is A2.22 = 
2.5 ± 0.5. Extrapolating this extinction to the optical is known to 
be highly uncertain, but we estimate Ay = 45 ± Mc81 is 
therefore one of the most heavily reddened clusters known, with an 
extinction comparable to that of the Galactic Centre. 

We estimate the distance to the cl uster from the radial veloc- 
ity of the surrounding molecular cloud. ICaswell & Haynes |l987t) 
studied the radio recombination line emission of the two clouds of 
ionized gas either side of the cluster, G338.4+0.2 and G338.4+0.1 
(see Fig.Q}, finding velocities relative to the local standard of rest of 
ui sr =-29km s _1 and -37 km s _1 respectively. In addition, a number 
of massive young stellar objects (YSOs) and compact Hll-regions 
have be en detected in the region by the Red MS X Source (RMS) 
Survey jHoare et al]|2005l : lUrquhart et"ai]|2007allbh . with an aver- 
age radial velocity of t>i sr = -35.1±2.8 km s . 

Co mparing this a verage «i 8I to the Galactic rotation 
curve of lBrand&Blitzl dl993h . using a Galacto-centric distance 
of 7.6±0.3kpc and a rotational velocity of 214±7 km s _1 



3 Using the value of a = —1.53 from Rieke & Lebo fskvl jl985l) . we find 
A2.22 = 3.5 ± 0.5, and A v = 30 ± 4 



dKothes & Doughertvll2007l, and references therein), we find near 
and far kinematic distances of 3.8kpc and llkpc. To resolve the 
near/far ambiguity, we refer to the Southern Galacti c Plane Survey 
(SGPS ) of Hi. These data were recently analysed bv lLemiere et al.l 
(2009), who found that for the two large Hll-regions G338.4+0.2 
and G338.4+0.1, Hi absorption components could be seen at sev- 
eral velocities up to the tangent point velocity of 130 km s _1 . This 
indicates the Hll-regions, and by association the YSOs and the cen- 
tral star cluster, lay beyond the tangent point. From this we con- 
clude that the clusters lay at the far-side distance of 1 1 kpc. At this 
distance, the uncertainty is dominated by the systematic uncertain- 
ties in the Galactic rotation curve, which by its nature is difficult 
to quantify. However, if we assume th at the system m ay have a pe- 
culiar velocity of up to ±20 km s~ _1 jRussei]||2003h . this gives an 
uncertainty on this distance of ±2 kpc. 

We can check this distance by calculating the absolute bright- 
ness of the WNL star and comparing to similar objects. Using 
the extinction calculated in the previous section and a distance of 
ll±2kpc, we find an absolute magnitude for Mc81-2 of M222 = 
— 7.1±0.6. By comparison , Galactic WNL stars are typically found 
to have M K = 5.9 ± 1.0 JCrowther etalfeOQfj) . Mc81-2 is there- 
fore somewhat luminous for its spectral type, though it is within the 
errors for other Galactic WNL stars. 



3.4 Cluster size 

In Fig. [6] we illustrate the physical extent of the cluster. The fig- 
ure shows the locations of all stars in the NICMOS field-of-view 
which are brighter than the 50% completeness limit (m 222 < 17), 
overlayed with those stars which have colours consistent with the 
cluster (i.e. mi6o — 771222 > 2.0). The figure shows that there is 
a clear overdensity of stars at the coordinate centre (defined as the 
position of star Mc81-2). 

To measure the size and morphology of this overdensity, we 
first made a map of the stellar density by dividing the field up into 
square bins of size 3" x 3" and counting the number of stars per bin. 
To reduce noise, this map was then smoothed, such that the effec- 
tive resolution of the map was 6" (illustrated in the bottom corner 
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Figure 6. An illustration of the spatial extent of the cluster. Red crosses 
denote all stars in the NICMOS field brighter than the 50% detection limit; 
while blue circles indicate stars with colours of mi60 — t*222 > 2.0. The 
contours indicate where the stellar density is 25%, 50% and 75% of the 
maximum. As the data were smoothed to make the contours, the circle in 
the bottom left shows the size of a resolution element. 



Table 3. Best fitting model atmosphere parameters for Mc81-2. 



Parameter 


Value 


T eff (K) 


36000 ± 1000 


T T =20 (K) 


38000 ± 000 


Voo (km s _1 ) 


1350 ± 100 





1.25 


A(H/He) 


0.75 ± 0.25 


log(L/L ) 


6.3 ± 0.4 


log(M/M eyr^ 1 ) 


-4.2 ±0.2 



of Fig. [§). This resolution size was chosen as a trade-off between 
spatial resolution and signal strength, though our results were ro- 
bust to changes in this parameter. The ambient stellar density was 
found by computing the background level of this map using the 
GSFC IDL routine SKY. Finally, we computed isodensity contours 
in this map at percentiles of the maximum stellar density. 

We defined the size of the cluster where the stellar density 
drops to 50% of its maximum value, which we deem to be roughly 
equivalant to its half-light radiufl Once deconvolved with the ef- 
fective spatial resolution, we find that the cluster has major and mi- 
nor axes of 29" x 18". At a distance of 1 lkpc, this corresponds to 
1.5 x 1.0 pc, and so is comparable to other young Galactic clusters 
which typical ly have diam eters between l-2pc (e.g. Trumpler 14, 
Westerlund l; lFigej2008h . 



4 Ideally, to measure the half-light radius one would measure the cumu- 
lative surface brightness out to a distance where it becomes asymptotic. 
However, the field-of-view of our observations is too small to do this. 
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Figure 7. Geneva evolutionary tracks from Meynet & Maeder ( 200(J). The 
x and y axes plot the H/He mass fraction and luminosity respectively, while 
the colour of the track indicates the effective temperature at that point in the 
evolution. The initial mass of each track is indicated at the top. Our derived 
parameters for Mc81-2 are indicated by the position and colour of the star 
symbol. Though the symbol intercepts the 120Mq track in (H/He)-L space, 
the temperature of the star means that the best fit is actually with the 6OM0 
track (see also Fig. [8). 
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Figure 8. Probability map derived from fitting stellar mass tracks to the 
luminosity, temperature and H/He of the WNL star Mc81-2. The age and 
initial mass of the best fitting model are indicated by the cross, and iso- 
probability contours are drawn at 67% and 50%. 



3.5 Quantitative spectral analysis 

To model the WNL star in Mercer81 and estimate its physical pa- 
rameters, we proceed as in lNaiarro et al. I fe004ll2009h . Briefly, we 
have used C MFGEN, the iterative, non -LTE line blanketing method 
presented by iHillier & Milled |T998) which solves the radiative 
transfer equation in the co-moving frame and in spherical geom- 
etry for the expanding atmospheres of early-type stars. The model 
is prescribed by the stellar radius, R, , the stellar luminosity, L* , the 
mass-loss rate, M, the velocity field, v(r) (defined by the terminal 
wind speed Voo and the wind acceleration parameter /3), the volume 
filling factor characterizing the clumping of the stellar wind, /( r), 
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and elemental abundances. Hilli er & Miiieil dl998l . [l999r) present a 
detailed discussion of the code. For the present analysis, we have 
assumed the atmosphere to be composed of H, He, C, N, O, Si, S, 
Fe and Ni. Observational constraints are provided by the H, K-band 
spectra of the stars and the dereddened F166W, F190N and F222M 
magnitudes. 

Given the extreme sensitivity of the H and AT-Band He I and 
Hell line profiles ratios in this parameter domain to changes in 
temperature, we estimate our errors in the temperature to be below 
1000 K. Likewise, the relative strengths between the H and He lines 
constrain the H/He ratio to be within 0.5 and 1.0 by mass. The 
error on _R» and hence on L, and M is dominated by those in the 
assumed distance and the slope of the extinction law. 

The best-fitting model is overplotted in Fig. [5] The model pro- 
vides a good fit to the features of the observed spectrum, with the 
exception of the BrlO line, which is blended with emission from 
NlV/ClV/OlV. This discrepancy is due primarily to the deficien- 
cies in the CNO IV model atoms, the correction of which is beyond 
the scope of the current work. The model's physical parameters are 
given in Table l3.4l and are typical for a late-type WN star. For com- 
pleteness, we list the temperature at an optical depth of r = 2/3, 
which is comparable to the star's effective temperature; and the 
temperature at r = 20 which is comparable to the hydrostatic tem- 
perature of stellar structure models. In the following Section we 
use these results to estimate the age of Mc81-3, and therefore of 
the cluster itself. 



3.6 Cluster age 

In order to contrain the age of the cluster, we make a quantita- 
tive comparison between the physical properties of the WNL star 
Mc81-3 and the predictions of stellar evolutionary models. Under 
the assumption that the star and the host cluster are coeval, we can 
then estimate the cluster's age. 

For this analysis, th e models we have chosen are those of 
IMevnet & Maederi §000) which are optimized for massive stars. 
In our method, we linearly interpolate these mass tracks at inter- 
vals of IMq and 10 5 yrs. For each point on each interpolated mass 
track, we then calculate the probability that there is a match be- 
tween the mass track and Mc81-3, based on the star's luminosity 
temperature T, and H/He ratio A(He) derived in the previous 
Section. For the star's temperature, we use the temperature at an 
optical depth of r = 20, since thi s is more comparab l e to th e hy- 
drostatic temperature calculated bv lMevnet & Maederi j2000l) . 

Figure [7] shows the inter-related behaviour of the three vari- 
ables L*, T and H/He for models with a range of initial stellar 
masses. We also plot the derived physical parameters of Mc81-3. 
The plot shows that, although the star's luminosity and H/He ratio 
place it on the 12OM0 mass track, the temperature of the star (illus- 
trated by the colour of the plotting symbol) more closely matches 
the 6OM track. 

We assume that the errors on L*, T and yl(He) are gaussian, 
and therefore the probability p of a match between Mc81-3 and the 
mass track of initial mass m at time t is given by, 

/ ,s TT ( (Mi-Oi) 2 \ 
P{m,t) = jjexp ^-^ — — 2 — '—j (2) 

where O is the observed quantity (either L t , T or yl(He)), 00 
is its associated uncertainty, and M is the corresponding quantity 
predicted by the model mass track. Each term is therefore weighted 
by its associated uncertainty. 



In Fig. [8] we plot how the probability varies across the 
2-D plane of mass and stellar age for the rotating models of 
IMevnet & Maedeil §000). The maximum probability (p = 0.75) is 
obtained for an initial mass of M* =62M© and an age of 3.7Myr. 
The morphology of the iso-probability contours are highly non- 
gaussian, so for the experimental uncertainty we cannot simply 
compute the standard deviation. Instead, we take the iso-probability 
contour at 50% and determine the minimum and maximum values 
of mass and age for that contour. In this way, for Mc81-3 we find 
M-,, = 62±?M and an age of 3.7±°'gMyr. Using the non-rotating 
versions of the same stellar evolution models produces a slightly 
different morphology to the probability distribution, with a reduced 
probability, but with a best-fitting age and mass that do not differ 
significantly from that derived using the rotating models. 

The value we obtain for Mc81-3's age can be understood 
through a simple qualitative analysis of the star's parameters. The 
high luminosity clearly favours high initial masses, and therefore 
a young age. In addition, the He enrichment indicates an object 
which is in an advanced evolutionary state, and so older than 
~2Myr, but younger than the total lifetime of a high-mass star, and 
so therefore younger than ~5Myr. 

3.6.1 The impact of binary evolution on our derived cluster age 

Throughout this analysis we have assumed that Mc81-3 has 
evolved as a single star. However, there is the pos s ibilit y that 
the star is i n an interacting binary: both ILandi et al] J2006l) and 
IFunk et alj J2007l) have detected X-ray emission from the centre 
of the cluster, which could be evidenc e of a colliding wind binary 
system. As shown by Eldridge (2009) for the case of 7 2 Vel, in- 
cluding the effects of binary evolution in the analysis can alter the 
derived age of a star. 

In the case of Mc81-3, the star's high luminosity places a 
strong constraint on the initial mass of the star, and hence on the up- 
per limit of its age. Though binary evolution can affect the surface 
abundances and temperature of a star, and prolong its lifetime, it is 
unlikely to increase the star's maximum luminosity by more than 
~0.1dex, which is governed pr imarily by the initial stellar mass (all 
other parameters being equal) (Eldridge et al. 2008). An exception 
to this would be if two stars merged to produce a completely rejeu- 
venated and more massive star. In the absence of any evidence for 
such an event in the history of Mc81-3, we maintain that the upper 
limit to the star's age is that derived in the previous paragraph. 

The lower limit to the cluster age may be reduced if Mc81-3 
is in an interacting binary system. Mass transfer from the primary 
to the secondary star may speed up the rate at which H is depleted 
from the primary's surface. In this case, we would underestimate 
the stellar (and hence cluster) age by using single star evolution 
models in our analysis. 

However, a simple morphological analysis of the nebula sur- 
rounding Mc81, and a comparison to similar systems, serves as a 
sanity check on our age estimate. The cluster is located at the centre 
of a cavity, which was presumably evacuated by the winds, ioniz- 
ing radiation and SNe explosions of the most massive stars in the 
cluster. At the periphery of the cavity evidence of further gener- 
ations of star formation is seen, which may or may not have been 
triggered by feedback from the cluster. This morphology is reminis- 
cent of other cluster + nebula systems such as G305, NGC 3603 and 
NGC 346 to name but a few. Ages of these other s ystems are com- 
monly found to be 2-4M yr dDavies et alj 1201 Q : iHaravama et alj 
l2008l : lBouretetal.ll2003h , and therefore are consistent with our es- 
timate for Mc81. 
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Figure 9. Top down view of the Galaxy, showing the locations of young star 
clusters. The colour of the plotting symbols indicate the extinction of each 
cluster, derived from their near-infrared colours and ass uming an extinc- 
tion la w slope of -2.0. The spiral arms are those defined by Cordes & Lazid 
120021) . Sites of multiple star clusters are indicated by a filled cicle sur- 
rounded by an open circle. 



Table 4. Young star clusters in the Galactic Plane. The visual extinction Ay 
has been determined homogeneously from each cluster's infrared colour 
excess E(H — K), using the relation A v ~ 19 X E(H — K), which 
follows from a NIR extinction law with slope a = —2.0. See the referenced 
papers for detailed error analysis on the distances. 



Cluster 


I (°) 


D (kpc) 


Av 


Ref. 


ril 806-90 


10 


8 7 


30 3 




V> J1J 


1 1 


4 5 


20 9 


2 




12 7 


4 7 




T. 

j 


M17 


1 5 


2 4 




4 


Mc9 


22.8 


4 2 


19 




W42 


95 d 


2 7 






RSGC1.. 5, Quartet 


26.0 


6.0 


19.0 


7,8,9,10,11,12 


W43 


30.8 


6.2 


39.8 


13 


CI 1900+ 14 


43.0 


12.5 


13.3 


14 


Mc20 


44.2 


3.5 


17.1 


12 


Mc23 


53.7 


6.5 


6.6 


15 


CygOB2 


80.2 


1.4 


1.9 


16 


h+x Per 


135.0 


2.3 


1.9 


17 


[DBS2003]45 


283.9 


4.5 


7.6 


18 


Westerlund 2 


284.2 


8.0 


7.6 


19 


Trumpler 14 


287.0 


2.5 


2.5 


20 


NGC3603 


291.6 


6.0 


4.7 


21 


Mc30 


298.8 


7.2 


10.5 


22 


Danks 1 & 2 


305.0 


4.0 


9.5 


23 


Mc81 


338.4 


11.0 


41.7 


This work. 


Westerlund 1 


339.5 


3.9 


9.5 


24 


NGC6231 


343.5 


1.8 


3.8 


25 


[DBS2003] 179 


347.6 


9.0 


15.2 


26, This work. 



3.7 Cluster mass 

For the cluster mass, it is difficult to make an accurate estimate 
without further spectroscopy of the stars in the cluster. We can 
however make a rough estimate of the cluster mass from the emis- 
sion line stars. If we assume that all the nine strong Pa emit- 
ters listed in Table [2] are WRs, then since the age we derive for 
Mc81 is roughly the same as that of Weste rlund 1 (Wdl) (3-5Myr, 
ICrowther et al.l2006l : lBrandner et al.l2008h which has 27 WRs, this 
suggests that Mc81 may be a factor of ~3 less massiv e than Wdl. 
As most estimates of Wdl's mass are around 1O 5 M0 dClark et all 
l2005t iBrandner et al.l r2008h this implies that the mass of Mc81 
is a few x 1O 4 M0 . We stress however that this is only an order- 
of-magnitude estimate; a more precise measurement of the cluster 
mass awaits further analysis of its stellar population. 



4 DISCUSSION 

4.1 Location in the Galaxy 

With the many recent discoveries of young star clusters in the 
Galaxy, we can now begin to build up a picture of the Galaxy's re- 
cent cluster formation. In Fig.|9]we plot the locations of all known 
young Galactic clusters with distances from the Sun greater than 
~2kpc. All clusters in the plot are thought to have masses in ex- 
cess of 1O 3 M0 and ages <j20Myr. The references for each data- 
point are listed in Table [4] We have colour-coded each data-point 
according to its visual extinction, which we have calculated in a ho- 
mogeneous way from each cluster's E(H — K), measured either 
from the references listed in|4]or from 2MASS photometry, and an 
extinction law slope of a — —2.0 (see Sect. 13. 3t . 

Figure [9] shows that there are now a significant number 
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of reddened star clusters known at Galactic longitudes between 
10> I >50. This line-of-sight corresponds to the tangent of the 
Scutum-Crux arm, as well as the near end of the Galactic Bar 
(Benjamin et al. 2005). Since one would expect the star forma- 
tion rate to be comparatively high in this location of the Galaxy, 
it is also reasonable to expect it to be rich in young star clus- 
ters. Indeed, the region hosts five known YMCs, plus a substan- 



episode of ~ 10 6 Mp) around 20Myr ago 


Garzonetal.Jl997 


Lopez-Corredoira et al. 1999; Figer et al. 2006 


; Davies et alj2007 


Clark et alj2009b; Neeueruela et alj2010l. 201 ll). 



However, less is known about the opposite side of the Galactic 
Centre and far end of the Bar. There are two possible reasons for 
this: firstly, the larger distance, high extinction and larger number 
of foreground stars (due to the intervening Galactic Bulge) make 
it more difficult to pick out clusters in by-eye searches. Indeed, it 
is unlikely that Mc81 would have been found were it not for the 
four bright foreground stars (see Sect. l3.lt . Secondly, the fact that 
no star clusters are known in this direction means that investigators 
are less likely to search this region. This is in contrast to the near 
end of the Bar, where the initial discovery of the cluster RSGC1 in 
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this region bv lFiger et all [2006) led to the subsequent discoveries 
of a further 4 clusters within the same complex^. 

From our distance estimate of Mc81, we can place this clus- 
ter close to where we suppose the far-end of the Bar may be, as- 
suming an azimufhal angle of 44° and a bar length of 4.4 kpc 
(Benjamin et al. 2005). Another cluster nearby, which may too 
trace the end of the Bar, is [DBS2003]179. The distance to this 
cluster is not well known, and is based on spectro-photometric 
distance estimates fo r stars with unknown luminosity classes 
jBorissova et al.l20 08). We have reassesed the distance to this clus- 
ter using a similar methodology that we have presented here for 
Mc81. Specifically, we assume that the cluster is physically asso- 
ciated to its nearby molecular cloud, and use the massive YSOs 
detected in the cloud to det ermine the systemic radial velo city. We 
then use the SGPS survey dMcClure-Griffiths et alj[200^ > to mea- 
sure the velocity spectrum of the intervening Hi gas to resolve 
the distance ambiguity. The average radial velocities of the YSOs 
(u lsr =-36±3 km s -1 ), combined with the Hi absorption which is 
seen up to tangent-point velocities of 130 km s~ , give a far-side 
kinematic distance of 9 kpc, again with a ±2 km s" 1 uncertainty 
to allow for deviations from the Galactic rotation curve (see Sect. 
33}. This is within th e errors of the spectr o-photometric distance 
of 7.9 kpc derived bv lBorissova et al] (2008). 

These two clusters - Mc81 and [DBS2003]179 - are then 
the first young star clusters to be discovered in this region of the 
Galaxy. In addition to these c lusters a group of giant Hll-regions, 
of which G338.4+0. 1 is one dRusseill 2003). suggest that this re- 
gion may be an active star-formation site, similar to the region of 
the Scutum-Crux tangent at the opposite end of the Bar. As such, 
future targeted surveys of this region may unearth a number of other 
such objects. 



4.2 Association with HESS 1640-465 

As was noted in the introduction, the Mc81 cluster is located only 
a few arcminutes from the TeV source HESS 1640-465, which is 
likely to be powered by a neutron star. If the two were associated, 
it would allow us to estimate the initial mass of the neutron star's 
progenitor. Here we discuss the possible association between the 
two objects. 

HESS 164 0-465 is known to have a counterpart source at 
GeV energies dSlaneetalJlaOloh. and in the X-ray, detected by 



uev energies (Maneetai. lulu), and m tne A-ray, d etected by 
Swift, XMM and C handra dLandi etail 120061: iFunk et al] 120071 ; 
iLemiere et al.l2009h . Analysis of the XMM X-ray spectrum yielded 
a column density of uh = 6.llQ g x 10 22 cm -2 or 3.6j^J'g x 
10 22 cm -2 , depending on whether a pow er-law or a bsorbed black- 
body m odel was used dFunk et ai1l2007l) . However, ILemiere et al] 
(2009) analysed the Chandra data and found that they required a 
higher column density of 1.4 x 10 23 cm~ 2 to fit the data. Assum- 
ing a standard calibration be tween tin and optical exti nction Av of 
n H = 1.8 x 10 2 Mvcm" 2 jPredehl & Schmittl 19951) , this implies 
an visual extinction of between 20 and 70 mags, depending on the 
model for the X-ray emission. Though the errors are large, this is 
consistent with our measurment of the extinction to the cluster of 
Av = 45 ± 15. 

From this evidence, it seems likely that HESS 1640-465 is as- 
sociated with the G338.4+0. 1 Hll-region surrounding Mc81. How- 



5 Though RSGC1 a nd RSGC2 were first re cognized as associations of stars 
by Bica et al. (2003) and Stephenson ( 1990) respectively, the nature of each 
cluster was not understood until later. 



ever, its connection with the cluster itself is not clear. If the cen- 
tral source of HESS 1640-465 is a neutron star (as seems likely), 
and the progenitor was born with the cluster but was ejected, there 
are two possibilities: either the progenitor was dynamically ejected 
from the cluster; or it received a kick from the supernova (SN). 
The location of HESS 1640-465 at the centre of SNR 338.3+0.0 
dGreenl 120041) provides circumstantial evidence against the lat- 
ter explanation, since this suggests that the progenitor exploded 
close to its present location. If the progenitor formed with the 
cluster but was ejected at a time t c j ago, the ejection velocity is 
~ 20(t c j/Myr) km s~ (assuming a projected distance of 22 pc). 
Therefore, it is entirely plausible that the progenitor star formed 
along with the rest of the stars in Mc81 and was dynamically 
ejected during the formation of the cluster. 

Finally, there is the possibility that HESS 1640-465 formed 
out of the same molecular cloud as Mc81, and at a similar time, 
but that the two formed independently of one another. The mor- 
phology of the G338.4 region suggests inside-out star-formation, 
with the ~3Myr old cluster in the centre and a series of YSOs 
and UC-HlI regions at the periphery of the surround ing cavity, 
which have ages of a few xlO yrs dDavies et al ] |2011bl) . The lo- 
cation of HESS 1640-465 does not fit this picture, since the pro- 
genitor star must have formed at least 2Myr ago, which is before 
the first SNe occured in Mc81. However, there are other known in- 
stances of 'multi-seeded' star formation, where collapse occurs at 
multiple causally- unrelated sites across the host GMC (e.g. W51, 
ICfark et alj2009al) . 

In summary, we conclude that HESS 1640-465 is likely as- 
sociated with the star-formation region of G338.4+0. 1. However, 
we are unable to make a definitive association with the star cluster 
Mc81, and so we are unable to use the age of the cluster to estimate 
the mass of the neutron star's progenitor, as we were in the cases of 
e.g. RSGC1 and CI 1900+14 dDavies et al.l2008ll2009h . 



5 SUMMARY 

We have provided a near-infrared photometric and spectroscopic 
investigation of the candidate star cluster Mercer 81 (Mc81). We 
find that that a highly extincted (Av = 45 ± 15) cluster exists in 
the field identified by Mercer et al. (2005), but that the bright four 
stars at the centre of the field are unrelated foreground objects. The 
cluster is located at the centre of a cavity in a large Hll-region in 
the direction of G338.4+0.1, with evidence of on-going star forma- 
tion in the periphery of the cloud. Our analysis of the cluster has re- 
vealed nine stars with strong Pa emission, one of which we identify 
spectroscopically as a late-type N-rich Wolf-Rayet star (WNL), in 
addition to a luminous early A-type supergiant. Via detailed mod- 
elling of the WNL star's spectrum we estimate an age for the cluster 
of 3.7l p'gMyr. Under the assumption that the stars with strong line 
emission are WRs, we have made an order-of-magnitude estimate 
of the cluster's mass of ^ 1O 4 M . 

From a kinematic analysis of the host cloud, we obtain a dis- 
tance to the host star-forming complex of 11 ±2 kpc. Our distance 
estimate therefore places the G338. 4+0.1 complex in the same re- 
gion of the Galaxy as the far end of the Galactic Bar. The recent 
detection of another star cluster close to this location, as well as 
other giant Hll-regions, suggest that this region of the Galaxy may 
be as active in star-formation as the opposite end of the Bar, where 
a ~ 10 6 Mq starburst event is known to have occurred in the last 
~20Myr. A targeted search of the far end of the Bar will likely un- 
cover many more young star clusters, though the high extinction, 
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large distance and dense stellar field of the intervening Galactic 
Bulge will mean that such a search will require considerable obser- 
vational effort. 
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